Statistically Deformable 2D/3D Registration for Estimating Post-operative Cup Orientation from a Single Standard AP X-ray Radiograph

The widely used procedure of estimating post-operative cup orientation based on a single standard AP X-ray radiograph is known inaccurate, largely due to the wide variability in individual pelvic orientation relative to X-ray plate. CT-based 2D/3D rigid image registration methods have been developed to measure post-operative cup orientation. Although encouraging results have been reported, their extensive usage in clinical routine is still limited. This may be explained by their requirement of having a CT study of the patient at some point during treatment, which is not available for vast majority of Total Hip Arthroplasty procedures performed nowadays. To address this limitation, this article proposes a statistically deformable 2D/3D registration approach for estimating post-operative cup orientation. No CT study of the patient is required any more. Compared to ground truths established from post-operative CT images, the cup orientations measured by the present technique in a cadaver experiment showed differences of 1.7±1.4° for anteversion and difference of 1.5±1.5° for inclination. When the present technique was evaluated on patients' datasets, differences of 2.2±1.3° and differences of 2.0±0.8° were found for the anteversion and the inclination, respectively. The experimental results, though still preliminary, demonstrated the efficacy of the present approac

Assessing the Accuracy Factors in the Determination of Postoperative Acetabular Cup Orientation Using Hybrid 2D-3D Registration

Single standard anteroposterior radiograph-based methods for measuring cup orientation following total hip arthroplasty (THA) are subject to substantial errors if the individual pelvic orientation with respect to X-ray plate is not taken into consideration. Previously, we proposed to use a hybrid 2D-3D registration scheme to determine the postoperative acetabular cup orientation and developed an object-oriented cross-program called "HipMatch.” However, its accuracy and robustness have not been fully investigated. To assess the potential factors that may affect the accuracy and robustness of the hybrid 2D-3D registration scheme in determining the postoperative acetabular cup orientation, a comprehensive validation study using a cadaver pelvis was performed. Nine X-ray radiographs taken from different pelvic positions relative to the X-ray plate and two computed tomography volumes of the pelvis with one acquired before the cup implantation and the other acquired after the cup implantation were used in the validation study. Potential factors that may affect the accuracy and robustness of the hybrid 2D-3D registration scheme were experimentally determined. Our experimental results demonstrate that (1) the plain radiograph-based method is not accurate; (2) the hybrid 2D-3D registration scheme helps to improve the estimation accuracy; (3) the hybrid 2D-3D registration scheme can robustly and accurately estimate the cup orientation even when a big portion of the radiograph is occluded; and (4) image resolution has minor effect on the estimation accuracy. The hybrid 2D-3D registration scheme is an accurate and robust method to measure exact cup orientation in THA. It holds the promise to be a valuable tool for clinical routine usage for providing evidence-based informatio

Fully Automatic Segmentation of Lumbar Vertebrae from CT Images using Cascaded 3D Fully Convolutional Networks

We present a method to address the challenging problem of segmentation of lumbar vertebrae from CT images acquired with varying fields of view. Our method is based on cascaded 3D Fully Convolutional Networks (FCNs) consisting of a localization FCN and a segmentation FCN. More specifically, in the first step we train a regression 3D FCN (we call it "LocalizationNet") to find the bounding box of the lumbar region. After that, a 3D U-net like FCN (we call it "SegmentationNet") is then developed, which after training, can perform a pixel-wise multi-class segmentation to map a cropped lumber region volumetric data to its volume-wise labels. Evaluated on publicly available datasets, our method achieved an average Dice coefficient of 95.77 $\pm$ 0.81% and an average symmetric surface distance of 0.37 $\pm$ 0.06 mm.Comment: 5 pages and 5 figure

Reconstruction of Patient-Specific 3D Bone Model from Biplanar X-Ray Images and Point Distribution Models

Reconstruction of patient-specific three-dimensional (3D) bone model from biplanar two-dimensional (2D) X-ray images and point distribution models (PDM) is discussed. We present a stable and accurate approach combining regularized morphing and shape deformation, and show its application to reconstruction of proximal femur. A novel image-to-model correspondence building method using directly the edge pixels detected from the 2D images and the apparent contour extracted from the 3D model is proposed to convert a 2D/3D reconstruction problem to a 3D/3D one, whose solutions are well studied. Quantitative and qualitative evaluation results on eleven cadaveric dry bones are given which indicate the validity of our approach

Computer assisted determination of acetabular cup orientation using 2D-3D image registration

Background: 2D-3D image-based registration methods have been developed to measure acetabular cup orientation after total hip arthroplasty (THA). These methods require registration of both the prosthesis and the CT images to 2D radiographs and compute implant position with respect to a reference. The application of these methods is limited in clinical practice due to two limitations: (1) the requirement of a computer-aided design (CAD) model of the prosthesis, which may be unavailable due to the proprietary concerns of the manufacturer, and (2) the requirement of either multiple radiographs or radiograph-specific calibration, usually unavailable for retrospective studies. In this paper, we propose a new method to address these limitations. Methods: A new formulation for determination of post-operative cup orientation, which couples a radiographic measurement with 2D-3D image matching, was developed. In our formulation, the radiographic measurement can be obtained with known methods so that the challenge lies in the 2D-3D image matching. To solve this problem, a hybrid 2D-3D registration scheme combining a landmark-to-ray 2D-3D alignment with a robust intensity-based 2D-3D registration was used. The hybrid 2D-3D registration scheme allows computing both the post-operative cup orientation with respect to an anatomical reference and the pelvic tilt and rotation with respect to the X-ray imaging table/plate. The method was validated using 2D adult cadaver hips. Results: Using the hybrid 2D-3D registration scheme, our method showed a mean accuracy of 1.0 ± 0.7° (range from 0.1°to2.0°) for inclination and 1.7 ± 1.2° (range from 0.0° to 3.9°) for anteversion, taking the measurements from post-operative CT images as ground truths. Conclusions: Our new solution formulation and the hybrid 2D-3D registration scheme facilitate estimation of post-operative cup orientation and measurement of pelvic tilt and rotatio

An integrated approach for reconstructing a surface model of the proximal femur from sparse input data and a multi-resolution point distribution model: an in vitro study

Background: Accurate reconstruction of a patient-specific surface model of the proximal femur from preoperatively or intraoperatively available sparse data plays an important role in planning and supporting various computer-assisted surgical procedures. Methods: In this paper, we present an integrated approach using a multi-resolution point distribution model (MR-PDM) to reconstruct a patient-specific surface model of the proximal femur from sparse input data, which may consist of sparse point data or a limited number of calibrated X-ray images. Depending on the modality of the input data, our approach chooses different PDMs. When 3D sparse points are used, which may be obtained intraoperatively via a pointer-based digitization or from a calibrated ultrasound, a fine level point distribution model (FL-PDM) is used in the reconstruction process. In contrast, when calibrated X-ray images are used, which may be obtained preoperatively or intraoperatively, a coarse level point distribution model (CL-PDM) will be used. Results: The present approach was verified on 31 femurs. Three different types of input data, i.e., sparse points, calibrated fluoroscopic images, and calibrated X-ray radiographs, were used in our experiments to reconstruct a surface model of the associated bone. Our experimental results demonstrate promising accuracy of the present approach. Conclusions: A multi-resolution point distribution model facilitate the reconstruction of a patient-specific surface model of the proximal femur from sparse input dat

Scaled, patient-specific 3D vertebral model reconstruction based on 2D lateral fluoroscopy

Backgrounds: Accurate three-dimensional (3D) models of lumbar vertebrae are required for image-based 3D kinematics analysis. MRI or CT datasets are frequently used to derive 3D models but have the disadvantages that they are expensive, time-consuming or involving ionizing radiation (e.g., CT acquisition). An alternative method using 2D lateral fluoroscopy was developed. Materials and methods: A technique was developed to reconstruct a scaled 3D lumbar vertebral model from a single two-dimensional (2D) lateral fluoroscopic image and a statistical shape model of the lumbar vertebrae. Four cadaveric lumbar spine segments and two statistical shape models were used for testing. Reconstruction accuracy was determined by comparison of the surface models reconstructed from the single lateral fluoroscopic images to the ground truth data from 3D CT segmentation. For each case, two different surface-based registration techniques were used to recover the unknown scale factor, and the rigid transformation between the reconstructed surface model and the ground truth model before the differences between the two discrete surface models were computed. Results: Successful reconstruction of scaled surface models was achieved for all test lumbar vertebrae based on single lateral fluoroscopic images. The mean reconstruction error was between 0.7 and 1.6mm. Conclusions: A scaled, patient-specific surface model of the lumbar vertebra from a single lateral fluoroscopic image can be synthesized using the present approach. This new method for patient-specific 3D modeling has potential applications in spine kinematics analysis, surgical planning, and navigatio

FACTS: Fully Automatic CT Segmentation of a Hip Joint

Extraction of surface models of a hip joint from CT data is a pre-requisite step for computer assisted diagnosis and planning (CADP) of periacetabular osteotomy (PAO). Most of existing CADP systems are based on manual segmentation, which is time-consuming and hard to achieve reproducible results. In this paper, we present a Fully Automatic CT Segmentation (FACTS) approach to simultaneously extract both pelvic and femoral models. Our approach works by combining fast random forest (RF) regression based landmark detection, multi-atlas based segmentation, with articulated statistical shape model (aSSM) based fitting. The two fundamental contributions of our approach are: (1) an improved fast Gaussian transform (IFGT) is used within the RF regression framework for a fast and accurate landmark detection, which then allows for a fully automatic initialization of the multi-atlas based segmentation; and (2) aSSM based fitting is used to preserve hip joint structure and to avoid penetration between the pelvic and femoral models. Taking manual segmentation as the ground truth, we evaluated the present approach on 30 hip CT images (60 hips) with a 6-fold cross validation. When the present approach was compared to manual segmentation, a mean segmentation accuracy of 0.40, 0.36, and 0.36 mm was found for the pelvis, the left proximal femur, and the right proximal femur, respectively. When the models derived from both segmentations were used to compute the PAO diagnosis parameters, a difference of 2.0 ± 1.5°, 2.1 ± 1.6°, and 3.5 ± 2.3% were found for anteversion, inclination, and acetabular coverage, respectively. The achieved accuracy is regarded as clinically accurate enough for our target applications